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Abstract. We consider a system described by a Fokker-Planck equation with a new type of momentum- 
dependent drift coefficient which asymptotically decreases as —1/p for a large momentum p. It is shown 
that the steady-state of this system is a re-generalized Gaussian distribution, which is a non-Gaussian 
distribution with a power-law tail. 

PACS. 05.20.Dd Kinetic theory - 05.20.-y Classical statistical mechanics - 05.90.+m Other topics in 
statistical physics 

1 Introduction distributions is the presence of an asymptotic power-law 

tail, Cai x ) ~ M - " -1 : ana - W q (x) ~ x 1 ^, respectively. 



Non-Gaussian probability distributions are frequently ob- 
served in a variety of systems such as physical, chemi- 
cal, economical and social systems. Known examples of 
non-Gaussian probability distributions are Levy cv-stable 
distributions, which can be defined by its Fourier trans- W K (x) oc ^— k(3x 2 + \J\ + K 2 f3 2 x 4 ^j , < 2) (3) 



There is another type of non-Gaussian distributions 
with asymptotic power-law tails, which is called a K-generalized 
Gaussian, 



formation as 



It has been originally studied in the context of statistical 



C%{x) = — [ dkexp[ikx - C\k\ a ], (0 < a < 2) (1) physics by Kaniadakis g]. This K-generalized Gaussian 
2vr J 

can be derived by maximizing Kaniadakis' K-entropy un- 



and Tsallis' g-generalized distributions PQ, 

W q (x) cc [1 - (1 - q)[3x 2 } ^ , (1< q < 3) (2) 



der appropriate constraints. This K-Gaussian reduces to 
the standard Gaussian, exp(— (3x 2 ), in the limit of k = 0. 
in the non-extensive statistical mechanics based on Tsallis' For a large value of x, the K-Gaussian obeys a power-law 



entropy 1 . A common key feature of the both probability as W K (x) ~ x 



2 Wada: A nonlinear drift 

The ^-generalized distributions have been shown to well 
explain, for example, the energy distributions of cosmic 
rays [2 , and the size distribution of personal incomes [3] . 
In a previous work [HIS], we have studied the asymptotic 
behavior of the K-generalized nonlinear Fokker-Planck(FP) 
equation, which steady-state is a ^-generalized Gaussian 
distribution. Furthermore a ^-generalized Gaussian is also 
derived [6] by generalizing the log-likelihood function in 
Gauss' law of error, which is an original method devel- 
oped by Gauss himself to derive a standard Gaussian. 

On the other hand, Lutz [7] has recently shown an 
analytic prediction that the stationary momentum dis- 
tributions of trapped atoms in an optical lattice are, in 
fact, Tsallis' g-generalized Gaussian Later, Gaeta [H] 
showed its invariance under the asymptotic Lie symme- 
tries. The prediction was experimentally verified by a Lon- 
don team [9] . This anomalous transport is described by a 
linear FP equation with a nonlinear drift coefficient, 

K°\P) = 5* (4) 

which represents a capture force with damping coefficient 
a, and this force acts only on slow particles whose mo- 
mentum is smaller than the capture momentum p c . 

A characteristic feature of this nonlinear drift is that: 
for a small momentum \p\ < p c , the drift is approximately 
linear K ol (p) ~ —p, i.e., it reduces to a familiar Ornstein- 
Uhlenbeck process; whereas for a large momentum \p\ > 
Pc it asymptotically decreases as K ol (p) ~ — 1/p. 
In contrast to most systems with power-law distributions 
which are often described by nonlinear kinetic equations 



iich leads to K-Gaussians 

|10j . the above process is described by an ordinary lin- 
ear FP equation. Consequently standard methods can be 
applied to the analysis of the problem. 

It is worth stressing that the Lutz analysis is not re- 
stricted to anomalous transport in an optical lattice, but 
can be applied to a wide class of systems described by a 
FP equation with a drift coefficient decaying asymptoti- 
cally as —1/p. 

In this contribution, we propose another momentum- 
dependent drift coefficient K(p) given by equation (|16[) . 
which also asymptotically decreases as —1/p for a large 
momentum \p\ > p c . We consider the process described 
by the linear FP equation with this drift coefficient K(p). 
Next section provides a brief review of ^-generalized ther- 
mostatistics and some properties of K-generalized Gaus- 
sian. In section three we consider an ordinary linear FP 
equation with the proposed momentum-dependent drift 
coefficient K (p) and a constant diffusion coefficient D. It 
is shown that the steady-state of the FP equation with this 
nonlinear drift coefficient K(p) is a ^-generalized Gaus- 
sian. The deformed parameter k can be expressed in terms 
of the microscopic parameters. In section four the asymp- 
totic behavior of the FP equation is studied. It is shown 
that the non-increase of the Lyapunov functional associ- 
ated with the FP equation. Then we numerically analyze 
the time evolutions of numerical solutions against different 
initial probability distributions, and show the asymptotic 
convergence of the numerical solutions to K-Gaussian. In 
section five we discuss the relation between f3 and the aver- 
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age energy in the parameter region that the mean-kinetic bility distribution w{p) 
energy diverges. The final section is summary. 



2 ^-generalized thermostatistics 

We first give the brief review of the generalized thermo- 
statistics based on K-entropy defined as 



S K = -h 



dp w(p) \n K w(p), 



(5) 



for a probability distribution w(p) of the momentum p. 
Here fee denotes the Boltzmann constant, and ln K (x) is 
the K-logarithmic function defined by 



ln K (x) 



x — x 
2k~ 



(6) 



The K-entropy S K is a real-parameter (k) extension of the 
standard Boltzmann-Gibbs-Shannon (BGS) entropy. The 
inverse function of ln K (x) is expressed as 



exp K (x) = exp 



1 



arcsinh(Kx) 



nx + VT+ k 2 x 2 ^J " , (7) 

and called K-exponential function. For a small value of 
x, the ^-exponential function is well approximated with 
exp(x), whereas a large positive value of x, it asymptot- 
ically obeys a power-law exp K (a;) ~ x x / K . In the limit of 
k — both lri K (x) and exp re (x) reduce to the standard 
logarithmic and exponential functions, respectively. Ac- 
cordingly the S K reduces to the BGS entropy. 

Maximizing the K-entropy S K under the constraints of 
the mean kinetic energy and the normalization of proba- 



S 

Sw ' 



(s K [w\ -P / dp y w (p) ~ 7 / dpw{p)J 



leads to a so-called K-generalized Gaussian, 



w ME (p) = a exp K 



(8) 



(9) 



Here 7 is a constant for the normalization, and depends on 
/3, which controls the variance of w ME (p). The parameter 

a and A are K-dependent constants, which are given by 

1 

'1 - 



a 



1 + K 



a = vT 



K , 



(10) 



respectively. 

The ^-generalization of free-energy was studied in [11] 
and given by 



F K 



(11) 



where 



dp ■ 



1 r 



(w ME (p))^+(w ME (p)) 



(12) 



The K-generalized free-energy F K satisfies the Legendre 
transformation structures, 



F K = U-^S K , ^{ ■ ' ("1 



where 



r°° n 2 
U = J dp^-w ME (p). 



(14) 



3 Proposed nonlinear drift coefficient 



Let us consider the linear FP equation 

~(K(p) W (p,t)) + L dp 



^w(p,t) = ~(K(p)w(p,t))+D^w(p,t), (15) 
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with a constant diffusion coefficient D and the momentum- 
dependent drift coefficient, 

ap 



we can express the steady-state as 



K(p) = - 



(16) 



1 



where a is a damping coefficient and p c denotes a capture 
momentum. Note that this proposed drift coefficient K(p) 
also behaves as — p for a small momentum \p\ < p Cl and 
asymptotically decreases as —1/p for a large momentum 
\p\ > p c as same as K ol (p) in anomalous diffusions in 
optical lattice [7]. We introduce the associated potential, 



p ( p 

V(p) = arcsinh ( — 



\Pc 



which is related with K(p) by 



K( P ) = -a— V(p). 



(17) 



(18) 



3.1 Steady-state 

Next we show the steady-state w s (p) of the FP equation 
with the nonlinear drift (|16p is a K-Gaussian. To this, the 
steady-state condition ■^w s (p) = in Eq. (IT5|) leads to 
d , . . K(p) d aV(p) 

T^ w ^ = nr = -d P -^- (19) 

In the last step, we used the relation lfT8|) . Substituting 
equation ()17|) and after integration we have 



\nw s (p) 



ctp 2, ( p 2 

- arcsinh — 

2D \pl 



const. 



then the steady-state becomes 

,-,2 



w s (p) oc exp 



a Pc ■ . I P 

arcsmn 



(20) 



(21) 



2D VPS 
By using the definition ([7]) of K-exponential function, and 
introducing the two parameters as 



2D a 

K= — 2' ^ = 7T' 
ape D 



(22) 



W S (P) 



ex Ph 



(23) 



where Z K is the normalization factor. We thus found that 
the steady-state of the FP equation with the nonlinear 
drift coefficient K(p) is nothing but a K-generalized Gaus- 
sian. 

Remarkably, the parameter k can be expressed in terms 
of the microscopic parameters, i.e., a,D,p c in the FP 
equation. This fact allows us to give a physical interpre- 
tation of the K-generalized distribution, as similar as q- 
generalized distribution in Lutz' analysis [7 . For example, 
in the limit of p c — > oo, the drift coefficient K(p) reduces 
to — p of the standard Ornstein-Uhlenbeck process, and 
the deformed parameter k of equation (|2"2"|) reduces to 0. 
This corresponds to the standard case in which the steady- 
state w s (p) is a standard Gaussian. 

Note also that the parameter j3 is expressed as the ratio 
of the friction coefficient a to the diffusion coefficient D, 
in analogy with the fluctuation-dissipation relation. We 
emphasize that the parameter f3 is not equal to an inverse 
temperature, because w s (p) is not an equilibrium state but 
a steady-state, for which temperature is not well defined. 

4 Asymptotic behavior 

We here study the asymptotic solutions of the FP equation 
with the nonlinear drift coefficient K(p). 
In a previous work [31[S] we have studied the nonlinear FP 
equation associated the K-generalized entropy, and shown 
the existence of the associated Lyapunov functional, which 
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characterizes a long time behavior of the process described 
by the FP equation. Similarly, the Lyapunov functional, 

C(t) = U[w] - — S[w], (24) 

a 

is monotonically non-increasing, i.e., 4;C K (i) < 0, for any 
time evolution of w(x, t) according to the linear FP equa- 
tion (jT5j) . In equation 



S[w] = — fee 
is BGS entropy, and 

U[w] = 



dp w(p, t) \nw{p, t), 



dp V(p)w(p,t). 



(25) 



(26) 



is the ensemble average of the potential V(p). 
The proof of the non-increase of C(t) is as follows: 



dC(t) 
dt 



dp 



dp 



d_ 

dw 



V(p) w H w In w 

a 

D , 



dw(p, t) 
dt 



V(p) H Qnw + 1) 

a 



dp yK^w + D-u, 
dp-\ -K (p) + D — In w } <0. (27) 

-oo a I " OP ) 

From the fist line to the second line we used the FP equa- 
tion (|15p . and in the last step we used the integration by 
part. 

Thus C(t) is non-increasing and consequently £(t) is min- 
imized for the steady-state w s (p) as 

min£(t) = lim C(t) = U[w s (p)} - — S[w„(p)]. (28) 
t— >oo a 

Note that the last expression is the free-energy associated 
with the steady-state w s (p) of the linear FP equation. In 
contrast to the K-generalized free-energy (fl3)) associated 
with the nonlinear FP equation [H[S], S is the standard 
BGS entropy and U is the ensemble average of the non- 
linear potential V(p) in the relation (|28l) . 
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Fig. 1. A typical time-evolution of w(p, i) from an initial prob- 
ability density with triangle shape. The microscopic parame- 
ters are set to p c = 1, D = 1, a = 4, so that k = 0.5, j3 = 4. 

4.1 Asymptotic convergence to / -Gaussian 

In order to study a long time behavior of the FP equation 
with the nonlinear drift coefficient K(p), we performed 
numerical simulations against different initial probability 
distributions. We used a variant of the numerical method 
originally developed by Gosse and Toscani [12] for the 
Cauchy problem an the evolution equation. For the details 
of the numerical scheme, please refer to [T2"] . 

A time-evolution of the numerical solution w(p,t) of 
the FP equation with K(p) is shown in figure 1. Note that 
the numerical solution w(p, t) seems to be asymptotically 
approaching to the K-Gaussian. In order to confirm this 
property, we fitted the numerical solution w(p, t) at each 
time t with the K-Gaussian 



a{t) exp K [-b(t)p 2 ] , 



(29) 



where a(t) and b(t) are fitting parameters. Then the time 
evolution of the function defined by 

w(p,t) 



a(i)exp K [-6(t>2 



(30) 
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are studied. If a numerical solution w(p, t) is perfectly fit- 
ted with equation (|29|) , the function r](p,t) becomes iden- 
tically zero. In figure 2 the time evolutions of rj(p,t) and 
w(p, t) against an initial probability distribution w(p, 0) 
with a triangle shape are plotted. It is obvious from this 
figure that the function r)(p, t) is decreasing to zero as time 
evolves. This fact shows that the numerical solutions are 
approaching to the K-Gaussian. 

5 The relation between f3 and the average 
energy 

The K-Gaussian of the steady-state w s (p) is not normal- 
izable for the parameter region of 2 < or equivalently, 
ap^. < D. As similarly as in the momentum distributions 
in an optical lattice, the physical meaning of this is that 
compared with the random momentum fluctuations D, the 
capture force ap\ is too weak to keep the particle around 
the bottom (p = 0) of the potential V(p). In other words, 
the potential is too shallow to capture a particle with the 
large random momentum fluctuations. 

Next let us turn our focus on the parameter region 
of 2/3 < \k\ < 2(D < ap 2 c < 3D), in which the second 
moment 

/oc 
dpp 2 w s (p), (31) 
-oo 

of the K-Gaussian becomes infinite. Consequently the mean 
kinetic energy, (p 2 ) /2m, diverges, and it is the hallmark 
of an anomalous diffusion. Lutz [7 showed that an explicit 
correspondence between ergodicity breaking in a system 
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Fig. 2. Long time behavior of the quantity r)(p,t) given in 
equation (|30|) and the numerical solutions w(p, t) against an 
initial probability distribution with triangle shape. The mi- 
croscopic parameters are set top c = 1,D — l,a = 2, hence 
k = l,p = 2. The number of calculated points are 101. The 
best fitted parameters a(t) and b(t) in (|29[) are indicated in 
each plot. Inset figure shows w(p,t) at each time t. 
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described by power-law tail distributions and the diver- where E is the energy of a system in thermal equilibrium 

gences of the moments of these distributions, i.e., Enscm- with the temperature T. However, as pointed out before, (3 

ble average and time average of the dynamical variable p n is not an inverse temperature since the steady-state w s (p) 

are not equal to each other in the infinite-time limit, when is, in general, a non-equilibrium state, in which the tem- 

the 2n-th moment of the stationary momentum distribu- perature is not well defined, 
tion for a system described by power-law tail distributions 
diverges. His analysis is also valid for the present study be- -» 

6 Summary 

cause both momentum dependent drift coefficient K ol (p) 

and K(p) have the same asymptotic behavior 1/p, and We have proposed a momentum-dependent drift coeffi- 



consequently both steady-states are non-Gaussian distri- cient which asymptotically decreases as -1/p for a large 
butions with power-law tails. momentum \p\ > p c . We have studied a system described 

Whereas the mean kinetic energy diverges in this way, b y the FP equation with this drift coefficient, and found 
in the same region, let us consider the following average that the steady-state is a ^-generalized probability distri- 
energy bution. We performed the several numerical simulations in 

order to study asymptotic behaviors of the numerical solu- 



d 



p—V{p))= 1 w s {p). (32) 

u-P I J-oc ( v\ tions against the different initial probability distributions, 



1 + 



Since 



l d / pax P exp K (-/3^) 



integrating by part, the r.h.s. of equation (|32|) becomes 



and found that these numerical solutions asymptotically 
approach to the K-Gaussian functions. 
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